INTRODUCTION 


Presented in this paper is a brief description of STRESS 
(Structural Engineering Systems Solver) which is a system for struc- 
tural analysis by digital computer. It consists of a language which 
describes the structural problem and a processor which produces the 
requested results. STRESS is a generd purpose system in the sense a 
that it is capable of RE ETESEEN Caenee structural types and 
situations. The input Lanugdge is problem-oridnted, i.e., the only 
problem description required lis in engtngering jrather than computer 
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language. : : ; | 
On the Recumptich that ¢erecdeh te not jineer, a 

few words congerning the gendral nature of the structural design problem 

appear to be ja order. For éxenpie, consider. simple building 

The ienbers of i tf cturali system may be 

of steel, reig : Axe. figidly 

connected at fhe jeintas mf objective qt design is to evolvd a 

structure-whigh will -suipport ithe. imposed jloads thout excess{ve stress | 

or deformatiog and with nec i dum economy. =. i 
The anaijeie.of the selkctoeny simple trade_in Fig. 1 raguires 

the determination of 63 distinct force and moment components. This is 

accomplished by the sdition of Ar both? tadder of eqaé#ions. Forty-two 


of these are classified as equilibrium equations. The remainder express 


frame shown 14 Fig. 
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the compatability of distortions between the various elements. The 
total set of equations may be subdivided such that analysis requires the 
solution of 21 simultaneous equations. It should be apparent that 
rigorous analysis of a more sizable structure (e.g., a 20-story building 
frame) requires an enormous amount of computation and data processing. 
The problem is further complicated by the fact that the deformation 
of the individual members and hence the compatability equations depend 
upon the size and elastic properties of those members. Hence design 
must be an iterative process each cycle of which involves a new analy- 


sis of the complete structure and a revision of the member sizes. 


Fig. 1 Rigid Building Frame 
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A typical STRESS problem description is shown in Fig. 3. Although 
the example is trivial it serves to demonstrate the simplicity of the 
STRESS input language. The input, which is completely shown in the 
figure, can be n jni*& matter of minutes. An analysis of this 
structure pe ie es would require approximately one hour. 
Thus, ev¢n for this very simple case, the use of a computer becomes 


economic -—~For! ranger structures the input 5 ded, only 
Need of Pay See 


by the ditional. nume data required and a fap 
The important point to be made in connection with ‘9 3 is that 
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the program a consiate,of, PogRNOSTIng terms such ag, Nexemese rei Bodog | ‘9d “member”, 


etc. An ae trained in structural analysis can learn the STRESS 
language in a few hours. He is then in a position to analyze by com- 
puter the majority of structures which he encounters in practice. In 
other words, the engineer who will make the design decisions is in 
direct communication with the machine on his own terms. By this means 
the use of computers in structural engineering becomes economical, not 
only for the large, complex problems as at present, but for the routine, 
day-to-day analysis which comprises the bulk of bnal 


network. Fenves joined Branin‘2) in formulating and sc the 
linear structural analysis problem using rere theory. 
analysis uséd dk grmikes is based on this wo BoA ack ibe th Bonde has 
more recently been derived more explicitly by ‘Guage: 

The electrical network, gonsists, of. buanches an nodes, with scalar 
quantities being associated with the properties and variables. The 
structural analogy to the branch is the member, the joint to the node. 
The unknown at a point in @ member is not a scalar, but a vector. The 


tniog redjons js avonday 1edmem edt .erodoev esa opts neidaizay iaiot 

2 yd tud ,colisanotane:s teenkl a yd gea ,bestaler al t0dmem el} at 

STRUCTURE PROIAEOPRQREM Oot0t ond aotatseniit } exugit .soltasiotace:) xit2em 
‘Tyee Ae ie ee oa diiw ,zed0 ed? 03 sedmem tdgterta = to bae ex0 mo1t 


_ MeMBER OF WERURHGtRLeaMeLD veri? odd To tetver eoxok Lexecea edT .aeowsed: ak 


myoaen aa a aes mogodiro us nk ainenoqgmoo eozet teeqli eerdt ‘to aselascoo ett 

METHOD STTPRDRGeI8ceS OT . exe ody twode patton wigeengnce toemos eerds baa 
ry sear — assfoo doad .8 x.) eade 1 obte 
Ue 


3 450, 2256 ec3 no etusnoqmoo edt Pps: 


§0 
0. 8 8 ak noid aersotanatd $e 


al ‘AK 105 Wo b0Mfe: ods eltdw ,eqy? ods to notsoaut a aed oh: aswoadan “to tedden 
oe : a bi te re meoe Yam Sivnex aid? tiguie ou betase sedW .20n et cotdsion 
SRE TDs UgeiOmn sbodten tnetettth aotsatuques kasd 20% tady al to0% edd tud 
ene Tomcene antes (anm .tiveet s ah .neqyy enutcrsse Jamiettib so2 bees 
e mneen 24 LOA REACTIONS GReeUACEMER Beets ods 49% nestiixw need evad amety 
SOLVE ati seqy edt So eno Uno gitssets done 
i eldaT ,etayians Yo bodiom seeattitn to saemecaigqath a yatrebtancd 
sutot 1¢6q atnenogees tofsev —— to toduns muaiaia od} awoda oats 
anoltiigp, ot egy? iaxitowita 1 .(1t) 
acoliaupe xie aso staylens edd tot 
ed} eexa stnetatenoo qaises ee. “0% beriupet tatot 10q 
ores ayawis exe I eldaT ai awoda ton a2 nenognoe Saemeneiqath bas so10% 
edd mort bettimo od Yam eoulay orem seed? sbereblenco ed ton been bas 
-anaizd edd mott beseleb ens Loo. bee aOR sakbnoqaers00 eds bas ‘wtodoev 
eagiq s 10% anolsieleb edt Uileods amedon aveds a8 out -anoltamot 
od Y1ae2el0n anoLzaupe avoous? Lumtn to redanse ed at viao son -omett 
Ro Jnebnegebat ei sat ZETq edt to ‘iis teomis aud ciao ab & ot bled evioa 


sexta Totoev jasmeseigqalh tabot edd th ow vino begales equ Lexusowrte 


_;~-formationys ten Ad toe 2 size 6 x 6. Kach column correspohits!'t6!1C° a Th abs 

" BY the e of ap pofient on .Skeyggght side of the equaticn, oat * s 
+ each is used for the computat 3 one of the components on the «@SS 4/¢ e 
nap ,ona a 


joint variables also are vectors. The member unknown at another point 
in the member is related, not by a linear transformation, but by a . 
matrix transformation. Figure 4 illustrates the force tranefomeataon SRUT DENTE 
from one end of a straight member to the other, with no forces 4 Peo naan 
in between. The general force vector for the three diemnsignal, stmaér 350 445561) 
ture consists of three linear force components in an orthdgosdl' Gi ten IO FE 


nome: {T @oVIGAUL: FW Samui 
and three nt components acting about the axes. The general, trans, ; ; OUT AM 


leftmajde. It can be shown that tigjdisplacement transformation is 5 eo ,Oene @ 
and equal to the transpose pf the inverse of the force tins i) 7 IGMOM 


Thy sjetwork conceytzellows us 4 readily deal conceptually with . e 
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number of unknowns is then a function of the type, sale tus cull ail 
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solution is not. When stated so simply this result may seem obvious: ; - cMA E 
eli KA : 


but the fact is that for hand computation different methods nai den | > 
used for different structure types. As a result, many Computer, PYO7; ; 4 | 
grams have been written for the linear analpiisdf 2iithed’ striktirés,?- = ane ss 
ar ae J9270 2ACTTIARH 232RFOF Ftajuaal 
each treating only one of the types listed Table 1. BVoge 
Considering a displacement or stiffness method of analysis, Table 1 
also shows the ainiaom number of unknown vector components per joint 
(SF). For structural typem anther tWAT ENE seeet FeBee the .qgplitions 
for the analysis can editiy be! Sb faets by bat Bie? the six equations 


per joint required for the space frame. By taking consistent axes the 
force and displacement components not shown in Table i are always zero 
and need not be considered. These zero values may be. omitted from the 
vectors and the corresponding rows and columns deleted fromthe trans- 
formations. Figure 5 shows schematically the deletions for a plane 
frame. Not only is the number of simultaneous equations necessary to 
solve held to a minimum, but almost all of the program is independent of 
structural type, related only by JF, the joint displacement vector size. 


Force Transformation for a Straight Member 


4. 


Fig. 


Reduction of Transformation for a Plane frame 
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Fig. 
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Table 1 Member and Joint Unknowns 
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| BOL! SLO iT 9nt0%_ . 
| STRESS ‘s “{ntended te be. an informative. ‘and easily usable struc- 


| tural design tool. The designer then must be able to specify his 
problem to the machine easily, rapidly and concisely. He should be 
able to specify the problem as he thinks of it, not in terms of how 

| . the machine solves it. He should be able gto peciiy a provies without 

x, performing any computations duniag data Preperation. . This [implies 

that the processor will deal a mech groye .4 nimprely high 

the genetation ‘and solutiba f the; qnalyags qustiogs.” ’ Kod qxémple® the ;* 

equations relate imbal ejiced Joint fences whigh ean die. 

variety of load types considered ‘by tthe dqsignen. fae 

operate on joint coordinates while the designer might refete geometry 

to bays and stories, or spans. in the process of generating and solving 


the equations, and in this pre-and post-processing the machine must 
deal with “a great ‘dilotnt df @atas ‘adetly in asbdy -f008. . She uktimber 

of arrays and their sizes are variable functions of the input data, 

the structural type and size. The form and features of the STRESS 
processor are related to these problems and a desire that the processor 
be a dynamic entity expandable by engineers. 
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Figure 6. System Block Flow Chart 
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Table 2 
Program Blocks 


PROCESS 


Translation 
Consistency check, Internal representation 
Compute member stiffness matrices 


Modify stiffness matrices for member end 
releases 


Process all types of raw load data into 
equivalent joint loads 


Rotate member stiffness matrices into global 
coordinates 


Generate symbolically structural stiffness 
matrix 


Modify stiffness matrix and joint loads for 
joint releases 


Solve, matrix equation for joint displacements 


Backsubstitute for other results and print. 
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Data 
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Figure 7. Core Memory: Layout 
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space for, reference and use the arrays without requiring the location 
of the array to be fixed or even constant during a part of the solutions 
process. Figure 7 shows the noemak oo re vehement for STRESS opera- 
tion. A very small area (about 300 words) of fixed location variables 
is included in upper core; Bart of this area contains codewords #iieh 
are used to reference the arrays. :Thes:cudbwords contain such 40-84/:5 
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maining: memovy down to: the: top: ef! progrempaénsists of a pool for 22.) 
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(rather than physical) input fields, performins,,dig tion look-up, 
binary conversion, etc. This routine is called for evety| logical data 
field during translation of input data, by translation rams written 
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‘Sample Problem Specification 


STRUCTURE .SAMPLE PROBLEM | 
NUMBER OF JOINTS 8 © 
NUMBER OF MEMBERS 8 


NUMBER OF SUPPORTS 3. epee at 


NUMBER OF LOADINGS 2 

TYPE -PLANE® FRAME™ 

METHOD. STIFFNESS... & 
“TABULATE FORCES, REACTIONS 
JOINT COORDINATES — «=: 3 


2X -240. “SUPPORT 
5X 0. § : 
8X 240. $ 

4h Y 2460. 

7X 260. Y 2460. 

3Y 420. | 

6 X 240. Y 420. 
MEMBER: INCIDENCES 

12 1.. 

25 4 

387 

oe 

5 4 7 

64 3 

776 

8 3 6 

“MEMBER: PROPERTIES 


8 PRISMATIC AX 10.. 12.300, - 
4 PRISMATIC AX 10. [2 300. 
MEMBER. PROPERTEES PRISMATIC 
1 AX 20. 12 200, | 

2 AX 20. 1Z 200. 

3° AX 20. $Z 200, - 

5.AX.4Q. 12 300. 

6 1Z 180. AX 2u. 

7 1Z°180. -AX 20. 

CONSTANTS E 30000. ALL 
LOADING 1 UNIFORM ALL BEAMS 
MEMBER LOADS 

8 FORCE Y UNIFORM -0.1 

4 FORCE Y UNIFORM -u.1 

5 FORCE -Y-UNEFORM =071 - ms 
LOADING 2 WIND FROM RIGHT | 
JOINT LOADS — - 
‘o FORCE X. -2U. 

7 FORCE X -20. 

SOLVE THIS PART 
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statements prior to this command aaa Acces cat uaa 
problem. SH ie 5 — 

For efficient use of time-sharing, the npat ‘ts typed ‘in aeing | 
the CTSS monitor input program in a form which @PRESS -can: acospt aud’. 
execute, The remote console is usd. for. gontroliing 't tie: ‘pyocessor - 
and for immediate correction of errors so as NOt Ato: delay: the design... 
Answers to the specified problem are show’ in. Tate 7 ae 

The results show the forces acting on the member ends. and: actin = 
on the joints. With the solution of the member end forces, ‘the member 
is statically determinate, so that the forces and deformations in the. 
interior of the member can be determined by elementary methods . up to 
now the development of STRESS has concentrated: om the: overall problem... 
We are now, however, attacking such problems as the interior forces to 
develop a more effective design aid. The joint loads on support: joints. 
represent the reactions. While the difference between the calculated 
joint loads and the applied joint loads gives a measure of oi solu- 
tion accuracy. ‘ a 

The engineer may then wish to alter the problem: for: has Hebelon= 
ing design. In most cases the alterations wilt be! a turdtion ef the - 
obtained results which were not known during creation of: the: Anput 
file. He might then describe the differences in “the: new problem to 
the processor and obtain results for immediate comparison: and. evalua- 
tion of the merits of the tact of the design, ‘Tabie: 5 stops ‘the 
changes necessary to analyze the same structure. _with,. new member proper- 
ties as suggested by the first analysis. Table 6 shows the effects of: 
the changes. 

The STRESS system is in a continuing ‘state of development. It is = 
expected that its capability will be extended: to. ine}ude eynamic ana- 
lysis, investigation of structural stability, and the behavior ‘of in- 
elastic structures. It is hoped that ultimately STRESS. wilh become 
part of a larger system which will be an aid to. automatic, structural 
optimization, 


Table 4 
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Sample Problem Results 


STRUCTURE SAMPLE PROBLEM 


LOADING 1 UNIFORM ALL 


MEMBER FORCES 


MEMBER JOINT AXTAL 

FORCE 
uf 2 10.545 
Z 1 -10.545 
2 5 38.982 
2 4 -38.982 
3 8 22.473 
3 7 -22.473 
4 a 1.229 
4 4 -1.229 
5 4 “1.846 
5 7 1.846 
6 4 12.161 
b 3 -12.161 
7 7 11.859 
7 6 -11.839 
38 3 2.594 
8 6 “2.594 


BEAMS 


SHEAR 

FORCE 
-1.229 
1.229 
O.481 
-U.481 
U.748 
-0.748 
1u.545 
13.455 
13.366 
10.634 
“2.594 
2.594 
2.594 
“2.594 
12.1b1 
11.339 


STRUCTURE SAMPLE PROBLEM 


LOADING 1 UNIFORM ALL BEAMS 
JOINT LOAUS 
JOINT X FORCE Y FORCE 
SUPPORT & 
2 1.2292 10.5447 
5 -U.4814 38.9819 
8 ~U.7478 22.4734 
APPLIED JO 
1 -U,QUU00 u.ugdd 
3 U.UU00 O.UuU0UU 
4 u.,UuuU -u.UUUU 
6 -U,00U0 U.v0UU 
7 U.QUU0U U.,U0UU 


STRUCTURE SAMPLE PROBLEM 
LUADING 2 WIND FROM RIGHT 


MEMBER FORCES 


MEMBER JOINT AXIAL 
FORCE 
11.195 
“11.195 
10.377 
“10.377 
@27..573 
21.573 
13.334 
713.334 
16.467 
-lo.4b7 
8.188 
“3.188 
-3.188 
8.188 
41.600 
~11.6u0 


OaAWNOOUMF SE WUWKRR Ee 
SCWONW EN EER NK ERD 


SHEAR 

FORCE 
713.334 
13.354 
“14.732 
14.752 
-11.934 
11.934 
11.195 
-11.195 
13.385 
713.385 
“11.599 
11.599 
“%.4U1 
8.401 
8.188 
-3.138 


STRUCTURE SAMPLE PROBLEM 
LOADING 2 W!SD FROM RIGHT 


JOINT LOADS 


JOUNT X FORCE Y FORCE 
SUPPORT R 
2 13.5343 11.1955 - 
5 14.7321 1u.3774 - 
te) 11.9339 “21,5727 - 
APPLIED JO 
1 u,UUUU -u.vUuu0 
3 U.UuUu0] u. 
4 -U.UUU2 v.UUUU 
o -2uU.0U0UL -U.UUUU 
7 -2U.0UU1 uv. 0000 


PART 1 OF PROBLEM COMPLETED. 


BENDING 
MOMENT 
-92.604 

~202.414 

44.045 
71.498 
65.663 
113.813 
202.414 

-551,690 
628.394 

-300.563 

“148.203 

-318.751 
186.750 
280.204 
318.751 

- 280.204 


MOMENT 
EACTIONS 
-92.o0044 
G4.0448 
65.6634 
INT LOADS 
-0, 0000 
U.,UQUU 
-U.uu0U 
-0, 0000 
“0.0000 


BENDING 
MOMENT 
-1776.267 
-1423.969 
-189U.313 
-1645,392 
-1b59.204 
-1l1l94.g41 
1423.969 
1262,902 
1434.174 
1778.188 
“1051. 684 
-1U36.215 
“583,247 
-928.867 
1036,215 

923,867 


MOMENT 
EACTIONS 
1776.2675 
1890.43127 
1ob9. 2038 
INT LOADS 
-0,QUU0 
-U.0UUU 
-uU.GOUU 
-u.QUUU 
-u.00U0 


Table 5 
Modification Specifications 


STRESS 1S READY FOR INPUT, 


modification of first part - second cycle for member sizes 


TYPE 
TYPE 
changes 
All i cl <r eee ae 


member properties 
TYPE 

1 iz 800.6 

TYPE 

2 iz 889.9 


TYPE wi brates Fe 


3 iz 800.6 
TYPE 
4h iz 583.3 
TYPE 
5 iz 800.6 


FYPE 


6 iz 446.3 

TYPE 

7 iz 339.2 

TYPE 

8 iz 446.3 

TYPE 

solve 

PROBLEM CORRECTLY 


ae aera ae PEs 


SPECIFIED. 


SOLUTION WILL PROCEED. 
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Table 6 
Modification Results 


STRUCTURE SAMPLE PROBLEM 


MODIFICATION OF FIRST PART - SECOND CYCLE FOR MEMBER SIZES 
LOADING 1 UNIFORM ALL BFAMS 


MEMBER FORCES 


MEMBER JOINT AXIAL SHEAR BENDING 
FORCE FORCE MOMENT 
1 2 10.982 =1,5729 -127,062 
1 1 -10.982 1.729 -287,964 
2 5 38,177 0.712 68.235 
2 4 -38.177 =0..712 102.705 
3 8 22.841 1.017 92.710 
3 7 -22,841 -1.017 151.376 
4 1 1.729 10,982 287.964 
4 4 71.729 13,018 -532.213 
5 4 -1,831 12,993 585.341 
5 7 1.831 11,007 -347,003 
6 4 12,166 ~2,848 -155.832 
6 3 “12,166 2.848 -356.873 
7 7 11,834 2.848 195.627 
7 6 711,834 72,848 317.079 
8 3 2,848 12,166 356.873 
8 6 “2.848 11.834 -317,079 


STRUCTURE SAMPLE PROBLEM 
MODEFLCATION OF FIRST PART - SECOND CYCLE FOR MEMBER SIZES 
LOADING 1 UNIFORM ALL BFAMS 


JOINT LOADS 


JOINT X FORCE Y FORCE MOMENT 
SUPPORT REACTIONS 
2 1.7293 10,9823 -127,0625 
5 -0.7123 38.1766 68,2353 
8 -1,0170 22.8411 92,7096 
APPLIED JOENT LOANS 
1 -0,0000 0,0000 0.0000 
3 0.0000 0.0000 0.0000 
4 -0,0000 -0,0000 -0, 
6 -0.0000 0.0000 -0.0000 
7 0,0000 0,0000 -0,0000 


STRUCTURE SAMPLE PROBLEM 
MODIFICATION OF FIRST PART = SECOND CYCLE FOR MEMBER SIZES 
LOADING 2 WIND FROM RIGHT 


MEMBER FORCES 


MEMBER JOINT AXIAL SHEAR BENDING 
FORCE FORCE MOMENT 
1 2 9.680 -12,.474% -1790,689 
1 ut -9.680 12.474 -1203.048 
2 5 11,910 “15,680 -2144.816 
2 4 =11,910 15.680 -1618.268 
3 8 -21.589 711.857 -1760.015 
3 7 21.589 11,847 -1083.169 
4 1 12.474 9.680 1203.047 
4 4 “12.474 ~9,680 1120,.039 
5 4 16.769 14,103 1590.743 
5 7 -16.769 -14,103 1793.879 
6 4 7.487 -11.384 -1092,514 
6 3 -7.487 11.384 -956,610 
7 7 -7.487 -8,616 -710,710 
7 6 7.487 8.616 -840.170 
8 3 11.384 7.487 956.610 
8 6 “11.384 -7.487 840.170 


STRUCTURE SAMPLE PROBLEM 
MODIFICATION OF FIRST PART ~ SECOND CYCLE FOR MEMBER SIZES 
LOADING 2 WIND FROM RIGHT 


JOINT LOADS 


JOINT X FORCE Y FORCE MOMENT 
SUPPORT REACTIONS 
2 12.4739 9.6795 -1790,6894 
5 15.6795 11.9096 -2144,8154 
8 11.8466 -21.5892 -1760,0153 
APPLIEN JOINT LOADS 
1 0.0000 0.0000 -0.0000 
3 0.0000 9.0000 0.0000 
4 Qa,000a0 -0,0000 0.0000 
6 -20.0000 -0.0000 9.0000 
7 -20.0000 0.0000 -0,0000 
PROBLEM COMPLETED, 
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